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Abstract. The problem of quantifying the difference between evolutions of an open 
quantum system (in particular, between the actual evolution of an open system and the 
ideal target operation on the corresponding closed system) is important in quantum 
control, especially in control of quantum information processing. Motivated by this 
problem, we develop a measure for evaluating the distance between unitary evolution 
operators of a composite quantum system that consists of a sub-system of interest 
(e.g., a quantum information processor) and environment. The main characteristic 
of this measure is the invariance with respect to the effect of the evolution operator 
on the environment, which follows from an equivalence relation that exists between 
unitary operators acting on the composite system, when the effect on only the sub- 
system of interest is considered. The invariance to the environment's transformation 
makes it possible to quantitatively compare the evolution of an open quantum system 
and its closed counterpart. The distance measure also determines the fidelity bounds 
of a general quantum channel (a completely positive and trace-preserving map acting 
on the sub-system of interest) with respect to a unitary target transformation. This 
measure is also independent of the initial state of the system and straightforward to 
numerically calculate. As an example, the measure is used in numerical simulations to 
evaluate fidelities of optimally controlled quantum gate operations (for one- and two- 
qubit systems), in the presence of a decohering environment. This example illustrates 
the utility of this measure for optimal control of quantum operations in the realistic 
case of open-system dynamics. 



PACS numbers: 02.30.Yy, 03.67.-a, 03.65.Yz 



Submitted to: New J. Phys. 



Environment-invariant measure of distance between evolutions 



2 



1. Introduction 

Control and optimization of open quantum systems is a rapidly growing area of 
theoretical [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15] and experimental 
[16, 17, 18, 19, 20, 21, 22, 23] research. This research is connected, to a considerable 
degree, to a large body of work that has been recently devoted to integrating concepts 
and methods of optimal quantum control [24, 25, 26, 27] into the fields of quantum 
information and quantum computation (QC) [28]. In particular, optimal control theory 
was used to design external control fields for generating quantum-gate operations in the 
presence of an environment [29, 30, 31, 32, 33, 34, 35, 36]. In this context, an important 
problem is to evaluate the distance between different quantum operations, especially 
between actual and ideal processes [28, 37, 38, 39, 40, 41, 42, 43]. One obstacle in 
the comparison of such processes is the reliance of many existing distance and fidelity 
measures on the initial state of the system or, to a lesser extent, on an integration 
or optimization over all possible states. Another difficulty arises when we consider 
the evolution of the quantum system of interest, in this quantum information 

processor (QIP), that is coupled to a (finite-dimensional) environment. The main 
obstacle is that the unitary time-evolution operator for the composite system (consisting 
of the QIP and environment) and the target unitary operator for the QIP have different 
dimensions. 

In this work, we develop a measure (originally proposed in [44]) that quantifies 
the distance between evolutions of a composite (bipartite) quantum system in such a 
way that only the effect of the evolution operator on one component (the sub-system 
of interest) is relevant in the measure, while the effect on the other component (the 
environment) can be arbitrary. Due to the invariance to the transformation on the 
environment, it is possible to evaluate the difference between unitary evolutions of the 
composite system (including the environment) and the sub-system of interest alone 
(excluding the environment), i.e., between unitary quantum operations of different 
dimensions. Therefore, the proposed measure enables us to compare unitary quantum 
operations (or, more generally, quantum channels [28]) corresponding to dynamics of 
ideal (i.e., closed) and real (i.e., open) systems. When the sub-system of interest is 
not coupled to an environment, the proposed distance measure reduces to a form that 
was previously used to evaluate the fidelity of closed-system unitary operations (see, 
e.g., [45, 46, 47]). Another important advantage of this distance measure is that it is 
independent of the initial state of the system. This property is important, since the 
fundamental objective of QC is to generate a specified target transformation for any 
initial state. 

The development of an appropriate objective functional is a crucial issue for control 
and optimization (for both classical and quantum systems) [48]. Since maxima or 
minima of the objective functional determine control optimality, the selection of the 
objective is just as important as that of the physical model and optimization algorithm. 
Therefore, constructing an objective functional that represents control goals in the best 
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possible way is one of the first tasks in any application of optimal control. In addition, 
the emerging studies of quantum control landscapes [49], for state-preparation control 
[50, 51], observable control [52, 53, 54], and unitary-transformation control [55, 56, 57], 
emphasize the profound effects that the choice of the objective functional has on the 
optimization [the control landscape is defined by the objective as a function of control 
variables, J' = J'{C)]. In particular, the relationship between the control landscape 
structure and optimization search efforts was recently analyzed [58, 59]. The measure 
we present here is well-suited for directly evaluating distances between general operations 
resulting from control of open quantum systems. The utility of this distance measure 
for optimal control of quantum gates in the presence of an environment is illustrated 
with a detailed numerical example. 

This article is organized as follows. Section 2 presents an equivalence relation 
between unitary evolution operators acting on a composite (bipartite) quantum system. 
This relation is based on the effect of these operators on the sub-system of interest, 
while the effect on the other sub-system (representing the environment) is not taken 
into account and thus can be arbitrary. The equivalence relation produces a quotient 
space that is a homogeneous manifold of the original space of unitary operators. In 
section 3, a quotient metric is defined on this space, resulting in a generalized expression 
of the distance measure. This measure can be used to quantitatively compare unitary 
evolutions of the entire composite system and one of its components (the sub-system 
of interest), e.g., unitary operations of different dimensions. Particular forms of the 
general distance measure can be obtained using different matrix norms. In section 4, 
the distance measure is computed explicitly using the Hilbert- Schmidt norm, producing 
an analytical expression that is calculated directly from the unitary operator of the 
composite system. In addition, several properties of this specific measure are discussed. 
Section 5 develops the distance measure using the induced two-norm (i.e., the maximum 
singular value). In section 6, we relate the distance measure developed in section 4 to a 
few typical fidelity measures defined for a general quantum channel (acting on the sub- 
system of interest) with respect to a target unitary transformation. It is shown that the 
distance measure and an average channel fidelity both depend on the trace norm of the 
same matrix (as the matrix norm increases, distance decreases and fidelity increases), 
and that the distance measure sets the lower and upper bounds of this fidelity. Section 7 
illustrates, using a numerical example, how the distance measure developed here can be 
directly incorporated into the optimal control of quantum operations in the presence 
of a decohering environment. In this application, the actual unitary evolution operator 
for the composite system (consisting of a QIP and its environment) is compared to 
the target unitary operator for the QIP alone, which is of a smaller dimension. This 
example demonstrates the utility of the distance measure in optimal control of unitary 
quantum gates for QC. Lastly, section 8 concludes the paper with a summary of the 
results and a brief discussion of future directions. For the continuity of the presentation, 
some mathematical details and most proofs are relegated to appendices. 
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2. Equivalence classes of unitary evolution operators 

Consider a finite-dimensional, closed, bipartite quantum system with a Hamiltonian 
(possibly time-dependent) of the form 

H:=H, + H, + H,^„ (1) 

where Hs and He represent the Hamiltonians for the system S (e.g., a QIP) and 
environment S, respectively, and i^int represents the system-environment interaction. 
Let Tis and Tie denote the respective Hilbert spaces of the system and environment, 
where Ug := dim{7is} and Ue '■= dim{7ic}- Thus, the Hilbert space of the composite 
system C is He '■= Tis ® He and n := dim{7ic} = '^s'^e- Let {\i)}, and {\i) ® 

be orthonormal bases that span the Hilbert spaces Hs, He, and He, respectively. 

For a given Hilbert space H, the set of admissible states, represented as density 
matrices on H, is denoted by V{H). As such, a density matrix p G V{H) is a positive 
operator of trace one, i.e., p > and Tr(p) = 1. Elements of V{H) correspond to either 
pure [Tr (p^) = 1] or mixed [Tr (p^) < 1] states. Density matrices for pure states can 
also be expressed as where represents a normalized vector in H (i.e., 

= 1). The group of all unitary operators on H is denoted by \J{H). 

Let U{t) G U{He) denote the unitary time-evolution operator of the composite 
system C, whose evolution is governed by the Schrodinger equation {h = 1): 

= -mmt), (2) 

with the initial condition U{t = 0) = Ic, the identity operator on He- Given an initial 
density matrix p{t = 0) G V{Hc), this state evolves in time as 

p(t) = f/(t)p(0)f/t(t). (3) 

For any state p G V(Hc), the state of the system S, with which system observables are 
calculated, is described by the reduced density matrix ps := Tre(p), where Tie denotes 
the partial trace over He- Thus, from (3), ps evolves as 

P3(t)=Tr,[f/(t)p(0)t/t(t)]. (4) 
If two composite unitary operators U,V & l]{He) are such that 

TTe{UpU^)=TTe{VpV^), VpGP(Hc), (5) 

then they will produce identical evolutions of the reduced density matrix p^. In this 
situation, U and V are said to be system- or S- equivalent, a relationship expressed as 
U ^ V- The condition described by (5) defines an equivalence relation [60] (in addition 
to being symmetric and refiexive, iS-equivalence is also transitive, i.e., if Ui ~ U2 and 
U2 ~ ^3, then Ui ~ t/3). The necessary and sufficient condition for iS-equivalence is 
stated in the following theorem. 
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Theorem 1. Let U,V e U(K)- Then Ti^ (UpW) = TT^{VpV^) for all density 
matrices p G ^^(Tic) if o-nd only if U = [Mg'?^ for some $ G U(7ic)- 



This theorem is related to a unitary invariance of Kraus maps [28, 61], which is 
described further in section 6. 

For any U G \J{TCc), let [U] denote the equivalence class of U, i.e., [U] = {(lLs®$)f/ : 
$ G U(7Yc)}- As a consequence of Theorem 1, every equivalence class is isomorphic to 
the closed subgroup Is ® U(^o) = {Is ® $ : $ G U(7^c)}, i-e., [U] is just a right- 
translation of Is CS> U(7ic) by U. Therefore, if only the evolution of S is relevant, then 
just the space of iS-equivalence classes, i.e., the quotient space Q := U(7ic)/ls ® U(He) 
needs to be considered, rather than the space of all unitary operators \J{Ttc)- While Q 
is not a Lie group like U(7ic), it is a smooth (right) homogeneous manifold of dimension 
(ng — l)nl [62]. Although not explored here, the fact that Q has the additional structure 
of a homogeneous manifold may be useful for further analysis (e.g., explorations of the 
landscape topology [57]) of the distance measure. 

3. Properties of metrics on the quotient space of equivalent operators 

A metric space consists of a pair of objects: a set Q and a real-valued metric or distance 
measure A. The ordered pair [fl, A) is a metric space if the following conditions are 
satisfied for any x,y,z E fl [60]: 

(i) Non- negativity: A(x, y) > 

(ii) Strict positivity: A{x,y) = if and only if x = y 

(iii) Symmetry: A{x,y) = A{y,x) 

(iv) Triangle inequality: A(x, z) < A(x, y) + A{y, z) 

In this work, f2 is the quotient space Q developed in section 2; elements of this space 
are the equivalence classes [U]. Let A be any metric on U(7ic), then the corresponding 
quotient metric A on Q is 



This expression can be used to obtain the quotient metric for various choices of A, for 
example, when A is the Hilbert-Schmidt norm distance, the induced two-norm distance, 
the geodesic distance with respect to some Riemannian metric, etc. [63, 64, 65]. 

The quotient metric (6) does not involve the partial trace over the environment, 
which is implicit in other distance (and fidelity) measures that use the operator-sum 
representation [28, 38, 39, 40, 41, 42, 43]. Moreover, the quotient metric is independent 
of any particular initial state, which is a very useful practical feature in applications 
to QC. In contrast to other measures of quantum gate fidelity, e.g., some of those in 
[28, 38, 39, 40, 41, 42, 43], A([[/], \V]) is directly evaluated by propagating the evolution 



Proof. See Appendix A. 



□ 



A{[Ul[V]) := min A((ls®$i)f/,(ls®<f2)l^). 



(6) 
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operator U{t) in (2), without making any assumptions about the initial state of the 
composite system C. This property allows one to consistently quantify the distance 
between unitary quantum operations (e.g., the actual and ideal operation) for any initial 
state. 

If A is a left-invariant metric on U(Hc) [i.e., A{WU,WV) = A{U,V) for any 
U,V,W e U(Hc)], then (6) becomes 

A([f/],[y]) =mmA(f/,(l,®$)\/). (7) 

In the context of a norm distance, (7) can also be expressed as 

A{[U],[V]) =Amin||f/-(l,®<l>)F||, (8) 

where A is a specified normalization factor and || ■ || is any left-invariant matrix norm 
on the space of n x n complex matrices, denoted as Ain (C). 

In the following sections, we introduce some useful properties of unitarily invariant 
quotient metrics on the homogeneous space Q. 

3.1. Closed-system target transformations 

A very typical situation in QC occurs when the target unitary transformation is specified 
for an ideal closed system (i.e., for S alone, excluding coupling to the environment). 
Let U denote the actual evolution operator for the composite system C and V denote 
the target quantum operation. If the target is defined for the closed system, then V is 
decoupled: V = Vs®Vc, where G U(7Ys) is the specified target unitary transformation 
for S and G U(?ic) is arbitrary. Because Vc is an arbitrary operation on it can be 
incorporated into $2 in (6): 

min A((l, ® $i)t/, (1^ ® ^2)% ® K)) = min AUI^ ® $i)[/, (1^ ® $2)(K ® %))■ 

Similarly, K can be incorporated into $ in (7) and (8). Therefore, for the closed-system 
target V^, expressions (6), (7), and (8) of the quotient metric are reformulated as 

A([f/], % ® IJ) = min A((l, ® $i)f/, {V, ® $2)), (9a) 

A([f/],[\/s®le]) =mmA(?7,(\/s®$)), (9b) 

and 

A([[/], % ® le]) = Amin ||[/ - (K ® <l>)||, (9c) 

respectively. From a practical point of view, these results mean that we can use 
the distance measure A([f/], [V"s ® Ue]) to quantitatively compare the actual evolution 
operator U G U(7ic) of the composite system C (of dimension n) and the target 
transformation V"s G U(7is) specified for the sub-system S (of a smaller dimension ng). 
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3.2. Chaining of unitary operations 

Consider a physical process composed of several components, such as the application 
of a sequence of n unitary operations, Un ■ ■ ■ U2U1, acting on the composite system, 
compared to the corresponding sequence of target unitary operations, Vn ■ ■ ■ V2V1, where 
Ui,Vi G U(7ic)- In this context, a metric or distance measure is said to possess the 
chaining criterion [12] if 



A{Um...U2UuV^...V2V,)<A{UuV,)+A{U2,V2) + ...+A{Um,Vrn). (10) 



Thus, the total error of the process is bounded from above by the sum of the individual 
errors A(f/j, Vj). This property is useful for analyzing the overall fidelity of a quantum 
algorithm, consisting of, e.g., a sequence of one- and two-qubit operations [28, 42]. 

Theorem 2. If A is a bi-invariant metric on U(7ic) fi-e., A{WiUW2,WiVW2) = 
A{U, V) for any U, V, Wi, W2 G \i{7ic)], then the chaining criterion (10) holds for any 



From this result, the chaining property holds for Vi G U(7-^s) ®U('^^e), i-e., the set of 
all decoupled or factorizable unitary operators in U(7Yc)- This condition is satisfied by 
any Vi that preserves the isolation of the system S from the environment £, a ubiquitous 
objective in QC. 

3. 3. Convexity 

Performing the minimization in (6) over $1, $2 ^ U('^^c) to compute A([f/], \V]) may be 
difficult and expensive in general. However, when the distance measure is determined 
by a left-invariant norm (7), the minimization can be recast as a convex optimization 
problem [66]. Relaxing the condition on $ from = Ic to < Ic produces 

such a problem since (a) any norm is a convex function, (b) the argument in the 
norm is affine in the optimization variable $, and (c) the set of matrices that satisfies 

< le is a convex set in $. To show this, note that the set of matrices that satisfies 

< Ue is equal to the set that satisfies ||$||2 < 1, where || ■ II2 is the induced two-norm 
described in section 5. Thus, ||6'$i + (1 - 6')$2||2 < 6'||$i||2 + (1 - 6')||"^'2||2 < 1 for all 
< 6* < 1, which follows from the triangle inequality. Expressing (7) [or (8)] as a convex 
optimization problem is desirable, provided the solution is on the boundary of < le 
(i.e., $ G U(7Ye)), because if a local minimum exists in a convex set, it is also a global 
minimum [66]. 

3.4. Stability 

Another important metric property is stability [-12]. A metric is said to be stable if 



f/i, U2,...,U^e U(K) and V,, V2,...,V^e U(7^s) ® U(K). 



Proof. See Appendix B. 



□ 



A{[U]®E,[V]^E) = A{[U],[V]), 



(11) 
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for any unitary operator S of any given dimension. Stability means that operations on 
an ancillary system that is not coupled to the original composite system do not affect 
the value of A, and hence A [12]. We do not prove stability for the general form of 
the distance measure in (6) or (8), but demonstrate it for the Hilbert-Schmidt norm 
distance in the following section. 



4. Computing the distance measure with the Hilbert-Schmidt norm 



In this section, the distance measure presented in section 3 is evaluated with the Hilbert- 
Schmidt norm (also referred to as the Frobenius norm), which is defined as [64] 



||M||hs := VTr(MtM) 



1 1/2 



i=l 



, \/MEMniC) 



(12) 



where cTj(M) is the ith singular value of M. The ith singular value, ai{M), is the square 
root of the zth eigenvalue of the positive square matrix M^M, where the singular values 
appear in descending order: crj(M) > (Ti+i(M) for all i. 

The general form of the Hilbert-Schmidt norm distance is, according to (8), 



AHs([f/],[V]) := A^min ||f/-(]l,®$)V^||Hs : $ e U(7^e 



(13) 

Using (12) in (13) and A„ = (2n)~^/^, the minimization over $ can be carried out 
analytically and, as shown below, we obtain: 

1/2 



AHs([f/],[V^]) 



n 



Tr 



where T G A^n,(C) is 

r := Trs {UV^) 
and II ■ IItt in (14) is the trace norm [64]: 



n 

I M II TV := Tr (^v^MtM) = ^(Ti(M), MM e Mr 

i=l 



(14) 
(15) 

(16) 



Note that in (15), Trs(-) denotes the partial trace over the system space Hg. 

Proof. The derivation of (14) begins with some algebraic rearrangement of (13), initially 
yielding 



AHs([f/],[^]) =min Jl- -Re<' Tr f/Vt (i, ® $t) 



n 



min \ 1 1 



1. 



n 



-Re< Tr, 



$t X Trs (UV^) 



min \ l 1 



-Re 



n 



Tr (r<l>t 



(17a) 

(17b) 
(17c) 
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Computing the singular value decomposition (SVD) of T yields T = WEX^ where 
W,X e U(Hc) and S = diag ((Ti, . . . , cr„J G .Mnc (I^)) with cTj > (Tj+1 > 0, for all i. 
Thus, 

Tr (r<l>^) = Tr (VTSXt^t) = Tr (SX^<l>^Vr) = Tr (SF) , (18) 
where Y = X^^W e U(7^c). Hence, 

Re [Tr (r$t)] = Re [Tr (SF)] = ^ a,Re (Yu) , (19) 

which, because Y is unitary, achieves the maximum of Tr(S) when Y = Ue- For a unitary 
matrix Y, as a result of orthonormality, = 1 for all i. Hence, < \yij\ < 1 for 

all i,j, and in particular < \yii\ < 1 for all i. Since S is diagonal and positive (with 
matrix elements ai), Re [Tr (SF)] = Y,i (^i^^ {Vu) < Zli (^■i\yii\ <Y.i'^i = Tr (S). Thus, 
the maximum of Re [Tr (SF)] is attained when Y = 1^. 

To show that Tr(E) = ||r||Tr, consider S = ly^TX = X''T''"H^, which implies that 
= X^^r^^rX. Because E is diagonal, the matrix is also diagonal, and X is the 
unitary matrix which diagonalizes the positive matrix F^^F. Therefore, it follows that 
= S = XtyftfX. Thus, 

max {Re [Tr (F$t)] } = Tr (S) = ||F||Tr, (20) 

and (14) results. □ 

From this proof it follows that 

mmjllf/- (1s®$)V||hs : = = mm | - (1, ® <I>)\/||hs : < Ic}, (21) 

and the minimum occurs on the boundary of the set $"1"$ < Ig, i.e., $ = WX'' G U(7ie), 
where X and W are the unitary matrices from the SVD of F. 

For the case where the target unitary transformation G U(7Ys) is specified for a 
closed system, results of section 3.1 can be used to obtain: 

AHs([f/], [V, ® le]) = f 1 - -lirik j , r = Tr, [U{Vj ® le)] . (22) 
4-1- Stability of the Hilbert- Schmidt norm 

Consider the stability [12] of the Hilbert-Schmidt norm, as expressed in (11): 

Ans{[U] ® H, [V] ®E)= A„„^ mm |||[f/ - (1, ® ^)V] ® H||hs : $ e U(K)}, (23) 

where S is any unitary operator and = dim{S}. For this norm, a tensor product of 
operators factors into the corresponding product of the norms of those operators: 

AHs([f/]®S,[V]®S)) = A„min{||[f/-(ls®$)V]||Hs} x ||S||hs/v^, (24a) 

= Xnmm{\\[U - (Is ® $)V]||hs}, (24b) 

= Ansm,[V]). (24c) 
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Therefore, with the appropriate normahzation factor, Ahs is stable. 

4-2. Unitary operators that are exact tensor products 

Suppose that U and V, the actual and target operators, respectively, are exact tensor 
products, i.e., 

U = U,®U^ and V = V^0 K, (25) 

where Us,Vs G U(7is) and Uc,Vc G U(7Ye)- Such a situation would occur if, e.g., there 
is no system-environment interaction (-ffint = 0) or the system and environment are 
decoupled at a final time t = tf. Then, from (15), 

r = Tr (U^^) U,V,\ (26) 

and 

lirilTr = |Tr (f/sK^) I Tr (^^VMIu^V,^^ = n,\TT (U^Vj) \ . 
The distance measure now becomes 



(27) 



AHs([f/s®f/c],[v;®K]) 



1 



l--|Tr(f/,y/)| 



1/2 



(2^ 



which is, as one would expect, completely independent of Uc and Ve. Indeed, as discussed 
in section 3.1, f/e and can be incorporated into $ in the quotient metric in (8) 
and thus simply replaced by Ig. Therefore, the matrix F in (26) can be replaced by 
r = Tr (UsVj^ lo [immediately leading to (27)], and the distance measure in (28) can 
be expressed as Ahs([^s ® le]) [K ® le])- 

The distance measure in (28) is related to other distance and fidelity measures used 
in previous works for optimal control of unitary operations [45, 46, 47] and analysis 
of the unitary control landscape [55, 56, 57] in closed quantum systems. For control 
landscapes, the analysis of critical points of AHs([f4 ® Ic], [Vg ® Ic]) with respect to the 
control field is especially important. 

Note that the distance in (28) is zero if and only if Us and differ only by a global 
phase, i.e., Ug = exp(i0)l^ for any real (p. Therefore, this distance can be considered as 
a phase-independent generalization of 

r 1 r 1^/^ 
Ug-Vg\\us=\l--Re[TT{UsVj)\j , (29) 



V2^ 

a frequently used distance [57] that is zero if and only if Us = Vs. 



5. Computing the distance measure with the maximum singular value 

Suppose a (possibly entangled) normalized pure state {ip) of the composite system C is 
acted upon by a unitary operator U. The output of this action, U\ip), can be compared 
to the output of the target action, V\ip) (modulo Ig ® The output state error e is 



e := [U-{ls®m] 



(30) 
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The maximum (over all normalized pure states) norm of the error is 

maxe^e = \\U - {ls®(!>)V\\l, (31) 

\i>) 

where || ■ II2 is the induced two-norm of a matrix [61], i.e., the maximum singular value 
of the matrix argument. | Now consider the distance measure based on the induced 
two-norm: 

A2([[/],[V^]) = A2mm|||f/-(ls®$)V^||2 | $ G U(K)}, (32) 

where A2 = 2^^. Although an explicit solution for A2([f/], [V]) is not presented in this 
work, we establish the following bounds: 

X2\\U - (Is ® < A2([f/], [V]) <X2\\U- (Is ® ^)V\\2, (33) 

where $ is the solution to the following optimization problem: 

$ :=argmin|||f/- (Is ® $)V"||2 : < Ic}- (34) 

Here, < Ue indicates that is a positive Hermitian matrix whose eigenvalues are 
at most 1. The optimization in (34) is convex for the same reasons as those discussed 
in section 3.3. However, there is no guarantee that the resulting optimizer $ is unitary. 
Since $ in (34) is less constrained than in (32), it follows that the lower bound in (33) 
applies. The operator $ is obtained from the SVD of $: 

^ = WSX^ '^=WXl (35) 

Thus, $ is a unitary approximation to $, and since it is not necessarily the optimal 
solution to (32), the upper bound in (33) follows. The lower and upper bounds in (33) 
will be close to each other if the singular values of $ are close to unity, i.e., if $ is close 
to a unitary matrix. 

For any matrix M G A1„,(C), the induced two-norm is bounded from above by the 
Hilbert-Schmidt norm [61]: ||M||2 < ||M||hs- Comparing the distance measure with 
these norms yields a similar relationship. Given U,V E U(7ic), let $hs be minimizer of 
AHs([f^], [V]). It follows that 

A2-iA2([t/], [V]) < \\U - (Is ® $HS) Vh 

< \\U - (Is ® $Hs) ^IIhs = K'^nsiPl [V]), (36) 

i.e., A2([f/], [V]) is bounded from above by {n/2y/^Ans{[U], [V]). 

I Equation (31) may also be expressed in terms of the operator norm [64]: 

(Is := supll [[/- (Is «) $)F] |V')ll2, where {ip\ip) = 1. In this context, || • II2 is 

W 

the standard vector two-norm. 
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6. Relating the distance measure to quantum channel fidelity 

Quantum operations, for example, a free or controlled time evolution (possibly in 
the presence of decoherence) or a transmission of information (possibly via a noisy 
connection), are generally described by quantum channels [28, 67]. In this section, we 
consider a common situation in QC, where the target quantum channel is a unitary 
transformation specified for a closed system (i.e., for S alone). We denote this target 
unitary transformation as \4 G U(7is) (cf. sections 3.1 and 4.2), and the Hilbert- 
Schmidt norm distance AHs([t^], ® Ic]) corresponding to the closed-system target 
is given by (22). We demonstrate that the distance measure AHs([f/], [Vg ® Ic]) and 
a typical measure of quantum channel fidelity both depend on the matrix norm of 
r = Trs [U{Vj ® le)] from (22). 

There are several ways to define the fidelity of a quantum channel with respect 
to the target unitary transformation [28, 37, 38, 39, 40, 41, 42, 43]. For example, the 
Uhlmann state fidelity, defined as [37, 38] 

^u(pi,P2) := IIVPiVPsII^ = {Tr [{y/plp2^/plY^^]^ , (37) 

where pi, p2 G ViTi), can be used to measure the effect of quantum channels on states. 
The Kolmogorov distance between two density matrices is [40] 

^k(Pl,P2) := ^IIPi - P2||Tr, (38) 

which bounds the Uhlmann state fidelity [40]: 

[l-Dk(pi,P2)]' <^u(pi,P2) < l-D^(pi,P2). (39) 

In this section, we develop a direct measure of fidelity for quantum channels based on 
T-a in (37), and establish bounds analogous to those in (39) with the distance measure 
Ahs- 

Let /C denote a completely positive and trace-preserving quantum channel [28], 
mapping the state ps to ps, where ps, Ps G V{Hs), with the following operator-sum 
representation (also known as the Kraus map [61]): 

Ps = /C[p,] ■.= J2K^PsKl (40a) 

i 

such that 

J2KjK, = l,, where K.eMnXQ- (40b) 

i 

The target quantum channel is the unitary transformation that maps ps to VsPsVj . 

We can compare the actual and target quantum channels, /C and V"s, respectively, 
using the Uhlmann fidelity of (37): 

} ■ (41) 



(\/sPsK^Ps) = {Tr (\/sv^V;^p,\/,VP^y/) 
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When the input state ps is pure, it can be expressed as ps = IV^s) ("^sh then (41) yields 

(v^mii^slKlPs) = {AlVjpsVM = 5^ K^slK^i^^lV^s)!' . (42) 

i 

To evaluate the proximity of the quantum channel /C to the target transformation V"s, it 
is convenient to use a measure that is independent of the input state. One such measure 
is the minimum pure-state fidelity ^p(/C, V^) [2s], which is obtained by minimizing the 
form (42) of the Uhlmann fidelity over the set of all pure input states: 

J-p(/C,Vs) :=min^,(l^,|^3)(V^s|K^Ps) =minV|(^3|K^i^^|^s)|'. (43) 

Calculating T^iJC^ K,) is not easy in general because the minimization in (43) is not 
convex. However, this problem can be related to a convex optimization, producing the 
following lower bound: 

Z(^,l^s) <-^p(/C,K), (44a) 

where 

T_{K,V,) :=minV|Tr (r/i^.p,)!'. (44b) 

I 

Both fidelities are contained in the interval [0, 1] and equal to one if and only if 
/C [ps] = KPsK;^ for all density matrices ps G V{T-Cs), i-e., if Ki = aiVs for all i and 
\aif = 1. For a given /C, ^{JC, Vs) is found via a convex optimization over all states 
Ps G P(7is), and hence, it can be efficiently obtained numerically. It is shown in [28] 
that T{}C, Vs) and Tp{]C, Vs) have the same pure-state minimum, but since the set of all 
density matrices is convex, while the set of all pure states is not, finding £(/C, V"s) is an 
easier optimization problem. 

If, instead of minimizing over all possible density matrices, ps in £(/C, Vs) is simply 
the so-called maximally-mixed system state pm,s := T^s/^s, this results in another variant 
of quantum channel fidelity, J-'m,s{J(^,Vs) [28, 41, 67]: 

^^,.(/c, v;) ■.= ^Y1 (^s^^O r • (45) 

This fidelity also evaluates the proximity of the quantum channel /C to the target 
transformation Vs [28] , but without state optimization. Interestingly, p^ s has the unique 
property of being the density matrix with the "shortest" maximum distance, under the 
Hilbert-Schmidt norm, to any other density matrix: 

Pm,s = argmin I max ||pi - P2||hs ) , Vpi,p2 G V{Hs), (46a) 

P2 V Pl / 
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( max IIpi - P2||hs ] = — (46b) 



where 



mm max \\pi 

This resuh suggests that pm,s is "centrally located" within the space of density matrices 
in the Hilbert-Schmidt geometry. 

Proof. See Appendix D. □ 

Relating U (the composite-system evolution operator) to /C (the quantum channel 
for the system S) requires the specification of an initial state of the composite system 
C. For simplicity, assume that this initial state is an uncorrelated tensor-product state: 

no 

p = ps® pc = ps®^C,u\^){v\, (47) 

u=l 

where, as before, p G Vijic), p > 0, and Tr (p) = 1. The reduced dynamics of the 
system S can be represented by the mapping in (40), where the Kraus operators K^^r 
are 



:= V^Tre [(Is ® W'){iy\)U] = (^8) 

i,i'=i 

with the unitary operator U expanded as 

Us Uq 

U=Y.H U,,>\^){i'\®\u){p'\. (49) 

•/I I T 

1,1—1 u,u' = l 

There exist infinitely many different sets of Kraus operators, {Kj}, that represent 
the same map K. (i.e., they evolve the system state ps in exactly the same way) [61]. This 
is related to the equivalence relation defined in section 2, namely that U and (Ig ® 
produce the same system evolution and hence, the same Kraus map. Moreover, any 
Kraus map for a fc-level system can be represented by a set of at most k"^ Kraus operators. 

The mapping from U and p to /C described above illustrates the dependence of the 
quantum channel K. (and hence, of the corresponding fidelities) on the initial state of the 
environment. Note that the exact treatment of the dynamics requires the propagation of 
the evolution operator for the composite system C, which is also required for computing 
A([t/], [V"s ® le]). Although approximations may simplify the calculation of /C, they 
may eliminate certain physical processes, e.g., the Markovian approximation limits the 
memory of the environment, preventing possible coherence revivals. However, even when 
no approximations are made, measures of fidelity for a quantum channel /C require (for 
the mapping from U to /C) the specification of the environment's initial state and thus 
are less general than the distance measure A([f/], [Vs ® le])- 

With (48), the fidelity J'ni,s(/C, V^) in (45) becomes 

1 .""^ 1 

^m,s(/C, K) = ^ C' l^uu'f = ^I|rv/Pl|||s, (50) 
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where pe ^ T^O^c) is the initial state of the environment, po = X]"=i C;^!^) ('^I- Also, we 
see that both the distance measure AHs([f^], ® Ic]) and quantum channel fidelity 
•^m,s(A^5 depend on the matrix F defined in (22). Specifically, increasing the norm of 
r will increase fidelity and decrease distance. As a concluding example, suppose that pe 
is the maximally mixed environment state, i.e., pm^e '■= ^e/^e- In this case, the quantum 
channel fidelity Tm,s{^, Vs) becomes 

^m,c(/c,r,) := — V |r,,f = — lirii^s, (51) 

further emphasizing the dependence of fidelity on the norm of F. 

We show that the inequality relation between the distance and fidelity of quantum 
states presented in (39) [40] also applies to the distance and fidelity of quantum 
operations given by AHs([f^], [V's ® Ij) and ^01,0(^5^)5 respectively. Specifically, we 
obtain the following lower and upper bounds of J-'^a^d^, Vs)'- 

- AnsiPl [Vs ® Ic])] ' < -^m,c(/C, V,)<1- AlsiPl [V, ® Ij). (52) 
Proof. See Appendix E. □ 



7. Applications: Optimal control of quantum gates 
7.1. Model open system 

In this section, we illustrate the use of the distance measure AHs([t/], [Vs ® Ic]) in the 
optimal control of one- and two-qubit gates relevant for QC. Specifically, this measure 
is used in numerical calculations to evaluate the distance between quantum operations 
generated in the framework of optimal control theory [25, 27] for a system coupled to a 
decohering environment [32, 33], and a target unitary transformation. We use a model of 
interacting two-level particles, which are divided into a QIP, composed of q qubits, and 
an environment, composed of e two-level particles. The qubits are directly coupled to a 
time-dependent external control field, while the environment is not directly controlled 
and is thereby managed only through its interaction with the qubits. The Hamiltonian 
for the composite system C, abbreviated a.s H = Hq + Hc(t) + H^^^,^ has the explicit form 

q+e q q+e q+e-1 

i=l 1=1 j=i+l i=l 

Here, Sj = (S'j^, Siy, Si^) is the spin operator for the ith particle (Sj = |crj, in terms of 
the Pauli matrices), Hq is the sum over the free Hamiltonians uJiSiz for all g + e particles 
[ui is the transition angular frequency for the ith particle), Hc{t) specifies the coupling 
between the q qubits and the time- dependent control field C{t) {fii are the corresponding 

§ Compare this expression to the abbreviated Hamiltonian in (1): Hq + Hc(t) ^ Hs + He- 



Environment-invariant measure of distance between evolutions 



16 



dipole moments), and H[^t represents the Heisenberg exchange interaction between the 
particles {'-fij is the couphng parameter for the ith and j'th particles). Hamiltonians of 
the form in (53) are often referred to as effective or spin Hamiltonians [68]. 

The evolution of the composite system C is calculated in an exact quantum- 
mechanical manner, by propagating the Schrodinger equation (2) for the Hamiltonian 
in (53), without either approximating the dynamics by a master equation or using a 
perturbative analysis based on a weak coupling assumption. This calculation produces 
the final-time evolution operator G U(7Yc) (an n x n unitary matrix), which is 
compared to the target gate operation Vg G U(7Ys) (an x Ug unitary matrix) via 
^Hs([t^tf], [Vs ® le]). Note that = 2^ = 2^ and n = 2^i+^l 

For our examples, we consider two different composite systems. In the first example, 
one qubit is coupled to a two-particle environment {q = 1 and e = 2), which can be 
modeled by a linear chain of particles with the qubit qi at the center, equally coupled 
to both environment particles 62 and 63: 

712 713 fr-A\ 

62 < — > qi < — > 63 , (54) 

where 712 = 713 = 7. Assuming only nearest-neighbor coupling in this configuration, 
we set 723 = 0. In the second example, two qubits are equally coupled to a one-particle 
environment {q = 2 and e = 1). This is modeled as a triangular cluster: 



63 

^13/ \T23 , (55) 

712 

qi < — > q2 

where the two qubits are denoted as qi and q2, the environment particle as 63, and 
7i3 = 723 = 7- Other system-environment configurations and corresponding optimal 
control results, using AHs([^tf], [Vs ® Ic]), are presented in [32, 33]. 

Parameters of the composite system are selected to ensure complex dynamics and 
strong decoherence: values of 7/0; are up to 0.02 and frequencies uji are close, but not 
equal, to enhance the system-environment interaction (see [32] for more details). In 
addition to setting h = 1, we introduce a natural system of units by setting /ij = 1 for 
all i and the qubit frequency ui = 1, implying that one period of free evolution is 27r.|| 



7.2. Optimal control algorithm 

The ultimate control goal is to decouple the system S from the environment S at a 
time t = t{ and simultaneously produce the target unitary operation Vs G U(7Ys) for 
S. Target quantum gates for the one- and two-qubit systems are the Hadamard gate 
Hg and controlled-not gate CNOT, respectively (both of these gates are elements of a 

|j For one-qubit coupled to a two-particle environment, the frequencies of the environment particles 
are: uj2 ~ 0.99841, 013 w 1.00159. For two qubits coupled to a one-particle environment, the frequencies 
of 52 and 63 are: 012 ~ 1.09159 and ui^ « 0.99841, respectively. 
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universal set of logical operations for QC [28]): 




and CNOT 



/ 1 \ 

10 

1 

10 



(56) 



A hybrid optimization method incorporating genetic and gradient algorithms [69] 
is employed to minimize the distance measure AHs([f^f], [^s ® le]) with respect to the 
control field C{t). After an initial optimization with a genetic algorithm, we remove 
the constraints imposed by the parameterized form of the control field to provide the 
potential for more flexible and effective control. An optimal control field is then found 
by minimizing the objective functional J'}^''\ given by 



J-f'(C) := AHs([f/*,(C)],[\/s®le]) 



a , 



CI 



(57) 



using a gradient algorithm (see [46, 32] for details of the gradient-based optimization). 
Here, denotes a particular closed subspace of L^([0,tf];M) (see Appendix F for 
details), and for the time-dependent Hamiltonian in (53), Utf : K^j U(7ic) denotes 
the map, defined implicitly through the Schrodinger equation (2), that takes a control 
field C(t) G to the unitary time-evolution operator Ut^ e U(7ic)- Thus, Kj^ is a 
Hilbert space of admissible controls, on which U{t;C) exists for all t G [0,tf] and all 



C e [70]. As such, 



[Vs] 



is the "dynamical" version of the distance 



measure Ahs, with an additional cost on the control field fiuence, where a > is the 
weight parameter for this cost. The gradient of J7^[^"^ is derived in Appendix F. 

7.3. Results 

Because < Ahs < 1 in general, it is convenient to define gate fidelities based on the 
respective lower and upper bounds of JFjn (.(/C, K) in (52): 



^nsiPt^l [Vs ® le]) := 1 - AHs([f/tJ, [Vs ® le]) 



(58a) 



and 



ns{[Ut,l [Vs ® le]) := 1 - A^s([f/,J, [Vs le]). 



(58b) 



Both fidelities are independent of the initial state and are evaluated directly from the 
evolution operator Ut^ of the composite system. These fidelities, computed for the one- 
qubit Hadamard gate, optimally controlled in the presence of a two-particle environment, 
are presented in figure 7.3 for various values of the qubit-environment coupling parameter 
7. For comparison, a fidehty based on the upper bound of the induced two-norm distance 
[as defined in (33)] is also presented: 



-^N ,, 1 2 



MPt,],[Vs0U) := [l-Asllf/i, -(K®$)||2] 



(59) 
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where $ is defined in (33)-(35). This fidehty was cast as a semidefinite program [66] 
and calculated from the final-time evolution operator Ut^ using the SDPT3 solver [71] 
and the YALMIP toolbox [72] in MATLAB. Since the evolution operators Ut^ produced 
via optimal control in this example are very close to the target equivalence class, the 
upper and lower bounds on the induced two-norm distance are nearly indistinguishable, 
although this is not true in general.^ 




y 



Figure 1. Gate fidelities (circles, also shown at a different scale in the inset), 
•^HS (squares), and J-2 (triangles) versus the coupling parameter 7, for the optimally 
controlled one-qubit Hadamard gate Hg. The qubit is equally coupled to a two-particle 
environment (g = 1, e = 2). Values of 7 range from to 0.02 in increments of 0.001. 

For the one-qubit system with no coupling to the environment (7 = 0), the optimal 
control field produces the Hadamard gate with the following values of the distance and 
fidelities: Ahs ~ 10"^, jFgg ^ 1 — 10"^^, jF^g ^ l — 2x 10~^. For the maximum coupling 
to the environment (7 = 0.02), we obtain for the optimally controlled Hadamard gate: 
Ahs ~ 0.0035, J^^fg ~ 1 - 10~^ ^hs ~ 0-993. For the two-qubit system with no 
coupling to the environment (7 = 0), the optimal control field produces the CNOT gate 
with the following values of the distance and fidelities: Ahs ~ 10~^, jF^g ^ 1 — 10~^, 
J^^s ~ 0.9998. For the maximum coupling to the environment (7 = 0.01), we obtain for 
the optimally controlled CNOT gate: Ahs ~ 0.02, J^^g ^ 0.9996, J^^g ^ 0.96. 

More significant than the actual fidelity values obtained in this example is the 
demonstrated ability of the distance measure AHs([^tf], [Vg ® le]) to quantitatively 
compare the evolution of a realistic open system to a target unitary transformation 

% For randomly selected unitary operators in U(7ic), the bounds in (33) are well separated and the 
minimizer $ defined in (34) is not unitary. 
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specified for an ideal closed system. Thus, the use of the distance AHs([f/tf], [Vs (g) Ic]) 
allows for the direct optimization of quantum operations in the presence of an 
environment, without specifying any part of the composite-system state. 

8. Conclusions 

We presented the novel and useful distance measure A([[/], [V]) for the quantitative 
comparison of unitary quantum operations acting on a composite quantum system, 
where the effect on only one sub-system (e.g., a QIP) is important, while the effect on 
the rest of the system (e.g., an environment) can be arbitrary. In practically important 
situations where the target operation is specified only for the sub-system of interest, the 
corresponding measure A([f/], [Vg® Ic]) can evaluate the distance between the evolution 
of the entire composite system and the target transformation (i.e., effectively compare 
unitary operations of different dimensions). This capability is especially desirable for 
measuring the distance between the actual controlled quantum process in a realistic 
open system and the target unitary transformation in its ideal closed counterpart. The 
measure developed in this paper does not require the specification of the initial state of 
any component of the system, is straightforward to calculate with the Hilbert-Schmidt 
norm, and can be used to define relevant state-independent fidelity measures. A fidelity 
measure independent of the initial state is extremely valuable for optimal control of 
quantum gates, which is a crucially important task in realistic QC. 

The quotient metric developed here is also applicable to other proposals for design 
and control of quantum operations. For example, it can be used for evaluating the 
fidelity of quantum gates in situations where multilevel encoding (MLE) of quantum 
logical states is employed (i.e., where each logical basis state of a qubit is encoded 
by multiple physical levels of a quantum system) [73] . Target unitary operations for a 
QIP equipped with MLE have a tensor-product structure similar to that in (25), and the 
logical equivalence can exist between different physical operations (i.e., different physical 
processes can result in the same quantum logical operation on the QIP). However, the 
environment is not explicitly included in the MLE formalism. 

Future research arising from the present work will involve the analysis of quantum 
control landscapes [55, 56, 57, 74, 75, 76, 77, 78] for objective functional based on various 
forms of the distance measure A([?7], [V^]), in particular, on the Hilbert-Schmidt norm 
distance Aiis{[Utf{C)], [Vg ® Ic])- It would be interesting to investigate the structure 
of these control landscapes and systematically search for robust control solutions that 
incorporate requirements on realistic physical resources. This search can be performed 
using a method called diffeomorphic modulation under observable-response-preserving 
homotopy (D-MORPH) [71, 75, 76, 77]. This method provides an accurate numerical 
tool for finding optimal controls. In particular, D-MORPH applied to closed quantum 
systems is able to identify optimal controls generating a target unitary transformation 
up to machine precision [77]. Extending the applicability of this approach to open 
quantum systems is an important goal, currently under development. Implementing 
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efficient numerical algorithms to calculate unitary time-evolution operators for larger 
composite systems is essential for practical optimal control simulations. Recent work has 
shown that such unitary maps can be efficiently constructed from appropriately designed 
state-to-state maps [79], utilizing the favorable state-control landscape structure [50, 55]. 
The possibility of expanding the method developed in [79] to the optimal control of open 
quantum systems is being explored. 



Acknowledgments 

This work was supported by the NSF and ARO. JD acknowledges support from the 
Program in Plasma Science and Technology at Princeton. Sandia is a multipurpose 
laboratory operated by Sandia Corporation, a Lockheed-Martin Company, for the 
United States Department of Energy under contract DE-AC04-94AL85000. 



Appendix A. Equivalence of unitary operations 

The theorem from section 2 is repeated here for convenience. 

Theorem 1. Let U,V e U(K). Then Ti^ (UpW) = TTc{VpV^) for all density 
matrices p G ViTi^) if and only if U = {Ms^ for some $ G U(7Ye)- 

Proof First, if t/ = (Is ® ^)V, then Tre (UpW) = Tre[(ls ® ^)VpV^ {t, ® <l>t) ] = 
Trc (ypV^^), so the "if" implication is proved. Next, to show the "only if" implication, 
suppose that Trc(f/pf/^) = Trc (V^pV^^) for all density matrices p. By complex 
linearity, Tre [UAU^) = Trc (VAV^) for all complex matrices A G A^„(C). Therefore, 
(5, Trc {UAW)) = (5, Trc {VAV^)) for all A G Mn{C) and all B G 7Wn,(C), where 

(fii,52) = Re[Tr(5|52)], Vfii, ^2 G A<„,,(C), (A.l) 
is the real Hilbert-Schmidt inner product on 7V1„^(C) [analogously for 7W„(C)]. If 

(5, Trc {UAU^)) = (fi,Trc (VAV^)), VA G Mn{C) and B G A<„,(C), (A.2a) 

then 

{Tt,*{B),UAU^) = {Tt,*{B),VAV^), ^AeMniC) and 5 G 7W„,(C), (A.2b) 

where Trc*(-) is the operator adjoint of Trc(-). Starting from (B, Ttc{A)) = {Ttc*{B), A), 
for all B G -A/1„,(C) and A G MniC), we show that Tre*(5) = B ^ 1^. 



ris n 



{B, Trc(A)) = E E i^^)^ ^^^^v (A.3a) 

i,k=l j=l 
Us n 

= EE(^'®le),«,.A,,,„ (A.3b) 

= (Trc* (5), A) Tre* (5) = (5 ® le). (A.3c) 
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As such, for all A e Mn{C) and B G Mn,{C), 

{B ®1,,UAU^) = {B ^1^,VAV^), (A.4a) 
{U^{B ® 1^)U, A) = {V^{B (g) 1,)V, A). (A.4b) 

Since (A.4b) is true for all A G 7W„(C) and B G A^„,(C), then 

U^{B ®1,)U = V\B ®1^)V, (A.4c) 
^ (fi®lc) (f/^^) = (t/W) (5(g)le). (A.4d) 

Finally, because (A.4d) is true for all B G MnAQ^ it holds for all W G U(Hs): 

(ly ® le) (f/V^^) = (f/^^) (W^® Ic), (A.4e) 

i.e., ?7y^ commutes with VT^Ue, for all W. As such, UV"^ is an element of the centralizer 
of \J(Hs) ® le [i.e., the subgroup of U(7ic) that commutes with U(7Ys) ® le, also formally 
defined in Appendix C], denoted as Cu(Hc)(U('^^s)®1c) [80, 81]. In Appendix C, we show 
that the centralizer of U(7is) ® Ic is Is ® U(7ic)- Thus, 

UV^ = ls®<^ els^UiUc). (A.5) 
So, we conclude that [/ = (l^ ® for some $ G V(Hc). □ 

Appendix B. Chaining of unitary operations 

After defining and establishing some mathematical objects and terms necessary for our 
work, we show that the chaining property from (10) holds for a unitarily bi- invariant 
quotient metric on the homogeneous space Q. 

Definition 1. Let G be a group and H be a subgroup of G. The normalizer of H in G, 
AfclH] = {(7 G G : gHg^^ = H}, is the largest subgroup of G such that H is a normal 
subgroup of Ng[H]. 

Lemma 1. The normalizer of U(7is) (g) 1^ in ^JiHc) is U(7Ys) ® U(7Yc) := (g) $ : 
\Ef G U(7Ys) and $ G U(7Yc)}- Likewise, the normalizer of Is ® U(7Yc) in U(7Yc) is 

U(Hs)®U(K). 

Definition 2. The centralizer of H in G, Cg[H] = {g G G : ghg^^ = h,yh E H}, is 
the subgroup of all elements in G that commute with all elements of H. Thus, Cg[H] is 
a subgroup of Mg[H] [80, 81]. 

Lemma 2. The centralizer of U(7^s) ® in U(K) is ® U(7^c) := { JZc ® $ : $ e 
U(7ic)}- Likewise, the centralizer of Is ® U(7Ye) in U(Hc) is U(7is) ® -Re- 
proofs of these lemmas are provided in Appendix C. The theorem from section 3 
is repeated here for convenience. 
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Theorem 2. If A is bi-invariant on U(7Yc); then 

A(f/^ . . . U2U,, V„,... V2Vi) < A {U^, V,) + A(f/2, ¥2) + ... + A{U„,, K,) (B.l) 

holds for any Ui, U2, . . . ,U.^ e U(Hc) and Vi, V2, ■ ■ ■ , e ^{U^) ® U(7^e)- 

Proof. We demonstrate that (B.l) holds for m = 2. This resuh can be generahzed to 
arbitrary values of m. If A is bi-invariant on \J(Hc), i.e., for any U, V, Wi, W2 G \J(Hc), 
A{WiUW2, W1VW2) = A(f/, V^), then, for any U2, 14, V2 G U(K), 

A{[U2U^l[V2V,]) =mm{A[U2U^,il,®<l>)V2V^]}, (B.2a) 

< mm {A[U2Uu f/2(ls ® $i)Fi] 

+ A [{U2ils ® (Is ® m2Vi] }, (B.2b) 

= A [t/i, (Us ® $ + mm I A [f/2, (Is ® $) ^2(13 ® $1)] } , (B.2c) 

< A[f/i,(ls®$i)Fi] +mm{A[t/2,(ls®$$2)^2]} 

+ A[(]ls®$2)V^2,V^2(ls®$I)]. (B.2d) 

The final result of (B.2) is obtained through repeated usage of the bi-invariance of A 
and the triangle inequality. Because this inequality holds for all $1, $2 G U(7ic), it holds 
for the $1 that minimizes A[Ui, (Is ® • Now, if V2 G J\fu(^Hc)[^s ® U(K)], i.e., V2 

is an element of the normalizer of Is ® U('^^e) in U(7Yc), then for each $1 G U(7Ye) there 
exists a $2 G U{He) such that A [(Is ® $2)V'2, V'2(ls ® '^'I)] = 0. In particular, such a 
$2 exists for the minimizing operator $1. Therefore, under these conditions, 

A([f/2f/i], [V2Vi]) < min {A[Ui, {1, ® } + mm {A[U2, (Is ® $$2)V"2] } , (B.3a) 

= A{[U^],[V,])+A{[U2],[V2]). (B.3b) 

□ 

From this result, the chaining property holds when V2 G Af\j(nc)[^s ® U(7ie)] = 
U(?^s) ® U('^^c), i-e., the set of all decoupled or factorizable unitary operators in U(7ic)- 
This condition is satisfied by any V2 that preserves the isolation of the system S from 
the environment S, a crucial target in QC. 

Appendix C. The normalizer and centralizer of U(7is) ® le in U(7ic) 

The lemmas from Appendix B are repeated here for convenience. 

Lemma 1. The normalizer o/U(7Ys) ® m U(7Yc) is U(7Ys) ® U(7Yc) := ® $ : ^ G 
U(Hs) and^ G U(K)}- T/ie normalizer of M^^V (He) U(K) «s a/so U(Hs) ®U(K)- 
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Proof. We determine the structure of the normahzer of U(?is) ® Ic in U(7ic); the 
structure of the normahzer of Ig C?) U(7ic) foUows from the symmetry of the tensor 
product. First, we define a fixed non-degenerate diagonal unitary matrix A G U(?is)- 
Note that if f/ G A/'u(>^^)[U(7^s) ® le], then in particular, t/(A(g) le)f^'^ e U(7^s) ® Ic- For 
simplicity, let f/(A ]le)t/^ = ^0 16, where V G U(7i;s)- Because U is unitary, V ® le 
and A ® Ue have the same eigenvalues (they are related by a similarity transformation), 
hence V and A must have the same eigenvalues, so there exists an G U(7is) such that 
QVn^ = A. Therefore, 

{n ® lc)t/(A ® lc)[(fi ® lc)f/]^ = A ® le, (C.l) 

demonstrating that {Q ® tc)U is an element of the stabilizer subgroup of A eg) le with 
respect to the group action of conjugation. Hence, 

^ V ' 

rts times 

i.e., {Q^lc)U is block diagonal with Ug blocks, each of which are rieXrie unitary matrices. 
So we have shown that any U in the normahzer can be written as ?7 = (fi ® tc)V for 
some n G U(7^s) and G U(K) © ■ ■ ■ © U(K). We can write 



V 



/ ^1 \ 



(C.3) 



where each Vi G U(7-^c) for all i. 

In order for U to be in the normahzer subgroup, it must be true that for each 
W G U(7^s), there exists an X G U(7^s) such that U{W ® le) = (X ® le)t^- If 
[/ = (r^t ® lc)V, we must have 1^(14^ (g) Ic) = (fiXfi^" © Ic)^^. This condition implies 
that WijVi = {QXQ)ijVj for all 1 < i,j < ris. Because W is arbitrary, for each i consider 
the case where Wn = 1 (and Wji = for j 7^ i), implying that is a scalar times Vi. 
Hence, V = diag [exp(i6'i), . . . , exp(i6'„J] © $ G U(7is) © U(?ie)- Therefore, we conclude 
that f/ = (fit ® 1^)V G ViHs) © U(K) and Aru(w,)[U(H,) © le] C ViH^) © U(K). 

For any f/ = ^ © $ G U(7^s) © U(7^e), f/(U(Hs) © lc)f/^ = ^U(Hs)^t ® 1^ = 
U(7^s) © Ic, so that 7Vu(7^,)[U(7^s) © Ic] D U(7^s) © U(K). Therefore, A/'u(w,)[U(7^s) © 
le] = UiK) © U(K) □ 

Lemma 2. T/ie centralizer of U(7is) © -^c U('^c) © U('^^e)- Likewise, the 

centralizer of 1^ © U(K) in U(7^c) «s U(7^s) © ^c- 

Proof. We consider the centralizer of U(7Ys) © Ig in U(7Yc), since the centralizer of 
le © U(7Yc) follows from the symmetry of the tensor product. Because the centralizer 
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must be a subset of the normalizer subgroup, each U G Cu(7^^)[U(7is) (g) Ic] must be of 
the form f/ = \1/ ® $ where ^ G U(7Ys) and $ G U(7Yc)- So, 

and ® ® lc)(* ® $)t = V ® 1^, VV G U(7^s)}• (C.4) 

Now, ® <I>)(V' ® le)(^ ® '^')^ = (^I/Vvl^^) ® Uc, so the only condition on \l/ and $ 
is that \l/ must he in the center of U(7Ys), i.e., \1/ = exp(i6')le and $ G U(7Ye)- Then, 
= exp{i9)l^®<^ = Ic® [exp(ie)$] G lc®U(Hc). Therefore, Cu(h,)[U(7^s) ® Ic] C 
le (8> U(7Yc)- Because every element of Ic ® U(7Ye) clearly commutes with U(7Ys) ® le, 
we conclude that Cu(We)[U(^s) ® Ij = Ic ® U(K)- □ 



Appendix D. Centrality of the maximally-mixed state 

We show that pm,s '■= ^s/^s is the density matrix with the "shortest" maximum distance, 
under the Hilbert-Schmidt norm, to all other density matrices, i.e., 

Pm,s = argmin ( max||pi - pallns ) ' Vpi, e 1^(7^8), (D.la) 

P2 \ pi J 

where 

/ ^ \ — 1 /T^ 1 \ 

min max ||pi — Pallns ~ • (D.lb) 

P2 \ pi J Us 

Proof. Let p2 = X]jCiK)(^l with < Ci < ■ ■ ■ < Cn < 1, and define the Hilbert-Schmidt 
distance as 

/(Pi) := \\p2 - PiIIhs = IIP2|Ihs + IIPiIIhs - 2Tr(p2Pi). (D.2) 
If pi = then /(pi) = ||p2||2jg + i _2Ci. Thus, max^, /(pi) > Hp^H^jg + i - 2Ci. 

Now, Tr(p2Pi) = XliO(^IPiK) for aUpi, where Y^MPiV) = Tr(pi) = 1, so Tr(p2Pi) > Ci- 
Also, 1 1 PiIIhs — general, which implies that /(pi) < ||P2||hs + ^ ^ ^Ci for all pi. 
Therefore, 

max ||p2 - Piling = IIP2|Ihs + 1 - 2Ci- (D.3) 

pi 

To minimize this expression over p2, we want to minimize ||p2||hs ^^"^ maximize Ci- 
Because Ci is the minimum eigenvalue of p2, its maximum value of l/ng occurs when 
P2 = ls/''^s- Likewise, ||p2||hs takes its minimum value l/rig at p2 = '^s/n^. So, indeed 
the argument of miup^ max^^ ||p2 — Pi||hs i^ Pm,s = Is/'^s and its value is (rig — l)/?7,s. □ 



Appendix E. Relating the distance measure to fidelity 

We demonstrate that the inequality relation between the distance and fidelity of 
quantum states presented in (39) [40] also applies to the distance and fidelity of quantum 
operations given by Ahs([^], [V^ ® le]) in (22) and J?-mc(/C, V^) in (51): 



Ahs([?7], [V, ® le]) < -F^,c(/C, V;) < 1 - Ahs([?7], [V; ® le])^ 



(E.l) 
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Proof. We first show that the bounds on T^ d^j Vs) follow from ||r||2 = maxj o"j < n^, 
where cxj is the ith singular value of F. Then we show that ||r||2 < holds. 

The upper bound holds if and only if ||r||^g < ns||r||Tr, or equivalently, — 
UsYlii^i- Assuming that maxj cxj < rzs, it immediately follows that Ylii^f — "^sYlii^ii 
which establishes the upper bound. 

With some algebraic rearrangement, the lower bound becomes ||r||Tr + 
(l/ras)||r||yg < 2y^||r||Hs- If ||r||2 < ^^s, which implies that ||r||yg < ns||r||Tr and 
||r||Tr + l/'^slirilHS — 2||r||Tr, the lower bound holds if ||r||Tr < v^H-'^Hhs- This is a 
known inequality for any matrix and thus establishing the lower bound. To prove this 
inequality, let a and r be two /c-dimensional vectors (where k < n^) with entries cTj, 
which are the singular values of F, and Tj = 1 for all z, respectively. It follows from 
the Cauchy-Schwarz inequality that ||r||Tr = Tli^i = ll*^!!! = ('^7 ^) — ||^||2||'r||2 = 
Hs < v^II-'^IIhs, where, for a vector, || ■ ||2 is the Euclidean norm (also 
referred to as the ^2 norm): 




11^112:= I >.^/ I . (E.2) 
To show that o",; < ng, consider 



= ||{Tr4f/(K^®]le)]}|^)||2, (E.3a) 
= llTrs [U{V^ ® le)(ls ® m] II2, (E.3b) 

where dim{|^/;)} = n^. and (V'l^) = 1- Applying the Cauchy-Schwarz inequality [64] to 
the argument of (E.3b) yields 

||r|^)||2 < v^llf/ [V^ ® Ic) (Is ® |V^))||2. (E.3c) 
Because the two- norm is bounded from above by the Hilbert-Schmidt norm [64], 

||r|^) II2 < v^ll {V} ® le) (Is ® m IIhs, (E.3d) 
= v^||1,®|V^)||hs, (E.3e) 
= n,{iP\iP)=n,. (E.3f) 

Since ||r|'(/^)||2 < is true for all normalized it follows that ||r||2 < "".s- n 
Appendix F. Optimal control objective functional: J'^^^ 

In section 7.2 we use the objective functional J7^[^°' in (57), which corresponds to a 
special situation (although ubiquitous in QC), where the target unitary transformation 
G U(7Ys) is specified for a closed system (i.e., for the system S). As discussed in 
section 3.1, in this situation the target unitary operator V G U(7Yc) should be replaced 
by the tensor product Vs ® le in the quotient metric, leading to the Hilbert-Schmidt 
norm distance AHs([t/], [K; ® Ic]) of (22). However, in this appendix, we consider a 
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more general case: an arbitrary target unitary operator V G U(7Yc), not necessarily the 
tensor-product one. The corresponding objective functional J'}^^ is 

J-f (C) := AHs([f/t,(C)], [V]) + ^\\C\\l^, (F.l) 

where the Hilbert-Schmidt norm distance AHs([f^], [V]) is defined by (14) and (15). 

Appendix F.l. Computing grad A([?7], [V^])hs 

Let H„ denote the space of nxn Hermitian matrices, i.e., H„ = {A G 7W„,(C) : A'^ = A}. 
We endow the linear spaces H„ and A^n(C) with the real Hilbert-Schmidt inner product 
in (A.l), and U(?ic) is given the Riemannian metric induced by this inner product. 
Define maps y, Z, and Ty by 



y:H„, y{Z):=y/l-n-^TiiZ), 

Z : Mn^C) ^ Z(r) := vW, (F.2) 

Fv^ : U(He) ^ MnSC) Tv{U) := Tr, [UV^ . 

For any fixed equivalence class [V^] G Q, let : U{Hc) — > M be defined by 

A|}3(t/) := AHs([t^], [V]). Then A|}3(f/) = yoZo Tvlu), and the differential of Aj^ at 
the point U G U(?ic) in the direction 6U G T[/U(7Yc) is given by the chain rule: 

duA^nk^U) = dzirymy o dryiu)Z o duTviSU)- (F.3) 

The differential and the gradient are related through the Riemannian metric by 
dc/A£3(5f/) = (gradA53(t/),5?7) for all 6U G Tt/U(K) [65]. Using (F.3), we can 
also write 

duA^nk^U) = (gradyjZ [TviU)] }, dr,iu)Z o duTviSU)), (F.4a) 
= ((d^ry)*o (dr,(t;)Z)*(grady{Z[ry(t/)]}), 6u), (F.4b) 

where gradyjZ [rv/(f/)] } denotes the gradient of y at the point Z [Ty{U)] and * denotes 
the operator adjoint. So, we can compute the desired gradient as 

gradAlg(f/) = {duryy[ {dr,.iu)Zy [gT8.dy{Z [TviU)] }]}. (F.5) 

To do this, we must find expressions for the gradient and two adjoint operators on the 
right hand side of (F.5). We develop each of these differentials separately. 



We begin by first computing giady. For y{Z) = yl — Tt{Z), the argument 
of the square root is an affine scalar function of Z, so dzy{SZ) = —TT{6Z)/[2ny{Z)]. 
The gradient is then given by dzy{5Z) = {grady{Z),6Z) = Re |Tr [grad y(Z)5Z] } = 
Tr [gTa,dy{Z)SZ~\. Hence, 
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We now turn to the problem of differentiating Z. Differentiating a square root in 
a commutative setting is straight-forward; we used this fact to differentiate y. But for 
matrices, non-commutivity may cause trouble. So, we begin by rewriting the definition 
of Z as = r^r, and differentiating both sides: 

ZdrZ{6T) + drZ{6T)Z = ST^T + T^T. (F.7) 

Define u : 7V1„^(C) to be the linear operator that stacks the columns of a square 

matrix to create a vector. If we let A^„,^(C) and C"'^ both be real Hilbert spaces with 
respective inner products {A,B) = Re [Tr (^t^)] and (x,y) = Re(xV), then z/ IS a 
linear isometry [60] between these two spaces: 

{A, 5) = 5^ Re (A^) Re (B,^) + Im (A,) Im (B,,) = {u{A), u{B)), (F.8) 

for all A,B & 7W„^(C). Hence u* = u^^ is the "matrix-ization" operator that returns a 
square matrix. 

From the discussion of matrix equations and the Kronecker product in [y2], we can 
rewrite (F.7) as 

[z^^i + i^z]u [drZ{6r)] = u{6r^r + r^^r) . (f.9) 

By construction, Z will be Hermitian and positive semi-definite; the transpose will 
also be Hermitian and will have the same eigenvalue decomposition as Z. Because of 
this, we know that Z"^ ® 1 + 1 ® Z will also be Hermitian and positive semi-definite (the 
eigenvalues of the Kronecker sum are just the set of all pairwise sums of the eigenvalues 
of Z). If we assume that Z is positive definite, then Z""" ® 1 -|- U ® Z will also be positive 
definite and this linear system will have a unique solution: 

drZ(5r) = u* \{z^ ® 1 + 1 ® z)"V(5r^r + r^^r)] (f.io) 

Let Cr : Mn^C) Mn^C) be given by Cr{A) = A^T + T^A. Then, considering 
Ain^^C) to be a real Hilbert space with the real Hilbert-Schmidt inner product, Cr is a 
linear operator, and we can find its adjoint: 

{A, C*rm = {CriA), B) = {AV + Ft^, B) , (F.lla) 
= Re[Tr(FU5 + A^F5)] , (F.llb) 

= Re|Tr {V B^)^ A + {V B)^ A ^ = {A,T {B + B^)) , (F.llc) 

for a given B E Mn^Q and all A e Mn,{C). Hence, C*^{B) = T{B + B^). Then, 

drZ = i^* o (Z^ O 1 + 1 ® Z) o z/ o £r, (F.12a) 

and 

drZ* = £* o I/* o (Z^ ® 1 + U ® Z)~^ o z/. (F.12b) 
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Using gra.dy{Z) = —l/[2ny{Z)] derived previously yields 

drZ*\gTady{Z)] = -—^—-C*^ou*o(Z^l + l^Zy'ou{l). (F.13) 

2ny[Z) ' 

Note that v* o {Z"^ ^l + l®^') "^ozy(l)is just the matrix solution to the problem 
ZA + AZ = 1, and therefore is A = \Z^^, since we are assuming Z is not just positive 
semi-definite, but positive definite. Then, 

d,Z-[grady(Z)] = (Z-) = -^TZ-\ (F.14) 

For the differential of Fy, we can extend Vy to Vy '■ A4n{C) 7W„^(C), by letting 
Tv{A) = Tts [AV^Y Then, Ty is a linear map, so the differential of Ty is just itself, 
i.e., dAt{6A) = ty{6A). For any 6U G Tc/U(K), duTy{6U) = dutv{SU) = ty{6U). 
We can find the adjoint of d^/Fy by 

{[duTyY{B),6U) = {B,duTy{6U)) = Re[TT,[B^duTy{6U)]y (F.15a) 

= Re|Tre[5tTrs (5f/V^)]}, (F.15b) 

= Re|Tr (lOfi"^)]}, (F.15c) 

= i ( [(1 (g) B)V - UV^ (1 B^) U] , 6U) , (F.15d) 

where the last step involves projecting (1 B)V onto Tc/U(7Yc)- Hence, (d[/Fy)* (B) = 
I [(1 ® B)V - f/^t (1 B^) U] for any 5 G 7Wne(C). 
Combining these components, we find that 

grad a12([/) = {duTy)* [ (dr^(c/)Z)* (gradyjZ [Fy(?7)] })] , (F.16a) 
= It, (R - UR^U) , (F.16b) 

where R = {1 ® {Ty{U)Z~^ [Ty{U)]})V = [l {nW^)]V and Ty{U) = nSW^ is 
the singular value decomposition of Ty{U). Since V G U(?ic) and fiVT^ G U(7Yc), 
we have R G U(?ic)- Note that R is the element of the equivalence class [V] that 
is closest to U under the Hilbert-Schmidt distance, i.e., R = (Is ® ^)V, where 
!> := argmin^. {||f/ - (1^ ® $)\^||hs}- 

Appendix F.2. Explicitly coupling the Schrddinger equation to the objective functional 

In this section, we consider the problem of coupling the Schrodinger equation to the 
kinematic cost function to arrive at a gradient fiow through the space of controls 
leading to an optimal control. To do this, we need to clarify the definition and 
geometry of K^^. Let rj : [0, tf] — >■ M, the "shape function" , be positive almost everywhere 
and continuous (and therefore bounded). Define the inner product on as 

(/.K:^fMd. (...) 
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where Kf^ is the set of (equivalence classes of) functions / G L^([0,tf];M) such that 
ll/IUtj < oo. Bilinearity, symmetry, and non-negativity of the inner product (■, ■)Kt^ are 
clear from the definition. When ?7(t) = 1, Kj^ is L^([0, tf];M) with the standard inner 
product. Other choices of 77 (t) change the standard geometry, moving undesirably- 
shaped functions far away from the origin, or out to infinity, where they are less likely 
to be the targets of an optimization on K^^. This point is discussed in more detail 
throughout this section. For the optimizations performed in this work, ?7(t) = sin(7rt/tf). 

To make this approach (and geometric interpretation) rigorous, we need to verify 
that Kjj is indeed a Hilbert space under the inner product (■, in (F.17). Suppose / is 
a measurable function such that ||/||Ktj < 00. Because rj is positive almost everywhere 
and bounded above, l/?7 is positive everywhere and bounded below by some (3 > 0. 
Then, 



IJ{^dt > m\i^ ■■= f^l' f{t)dt. (F.18) 



Hence, ||/||l2 < 00, so that / G L^([0, tf]; M), where || ■ \\l2 denotes the norm on [0,t{]. 
If f{t) 7^ on a set of non-zero measure, then /^(t) > on a set of non-zero measure 
[analogously for r]{t) and p{t)/ri{t)]. Hence, ||/||Ktj > 0, so we have non-degeneracy, 
and K^j is thus an inner product space. 

Suppose that {fk} is a Cauchy sequence in Kt^ [60]. Then gk := fk/^/v is a Cauchy 
sequence in L^, i.e., Hf^fe — S'^llLa = — /^HiKfj 0. Since L^([0, tf]; M) is complete, this 
implies that Qk ^ 9 & -Z^^([0, tf]; M). Let g := f/,yf], thus, 

ftf 



and 



/ ' f{t)/7]{t) dt= r g\t) dt = \\g\\l, < 00, (F.19a) 
Jo Jo 



II/.- /Ill, = r^^^^f^^^' ^F.i9b) 







/ ' [9k{t) - g{t)Ydt = \\gk - g\\l. ^ 0, (F.19c) 
^0 



so that / G IKjj and ^ / in the topology induced by (■, ■)Ktj- Hence, Kt^. with (-, ■)]^^^ 
is a complete inner product space, i.e., a Hilbert space. We have not simply effected 
a change in geometry, but we also have restricted the domain to a linear subspace of 
L2([0, tf]; R). If C G L2([0, tf]; M)/Ktj, then j}^ = 00, so that C is not a viable solution. 
This new geometry puts desirable functions near the origin and undesirable functions 
out near infinity, which makes this optimization different from the unshaped version. 

Now, consider the mapping Utf : U(7ic) introduced in section 7.2. The 

differential of the final-time evolution operator Ut^ with respect to the control field 
depends on evolution operators at all times in [0, tf]. It is shown in [77] that 

dcUt,{6C) = iUt,iC) r U^it; C)iiU{t- C)5C{t) dt, (F.20) 
^0 
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where /i is the dipole moment operator [e.g., /i = '^■fj.iSix in (53)]. With the usual 
Hilbert- Schmidt inner product as the Riemannian metric on U(7ic), we can find the 
adjoint of this differentiah for any A G T[/j^(f7)U(7ic) and all 6C G Tc^tf — ^tf, 



Re<! Tr 



A^Ut,{C) / U\t;C)fiU{t;C)5C{t)dt 







Re<^ Tr 



iA^Ut,{C)U\t-C)fiU{t;C) 6C{t)dt 



Rej Tr [iA^Ut,iC)U^ C)^iU{-- C)\\,5c) 



So that 



{[dcUt^rm (t) = -r/(t)Im{Tr [A^Ut,{C)U^{t;C)fiU{t;C)] }. 



(F.21a) 
(F.21b) 

(F.21c) 
(F.21d) 

(F.22) 



With jPiC) in (F.l), using previous arguments regarding the chain rule on 



gradients, we obtain 

(gradJ-f(C))(t) = \{dcU,,y{gT^dA^^^[Ut,iC)])]{t) + aCit), 



(F.23a) 



4nAl2[[/,,(C)] 
+ aC{t), 



Im( Tr I [R^ - U^tf, C)RU^{t{; C)]Ut,{C)U\t- C)fiU{t; C) 



lm(Tr{[i?tf/tj(C) - U^{tuC)R]U^{t; C)fiU{t;C) 



(F.23b) 

aC{t). 
(F.23c) 



Observe that grad Ahs is continuous on UiHc), and therefore bounded. It is shown 
in [77] that [dcUtf]* is uniformly bounded for all C G Kf^. In fact, ||(dct/tf)*|| < 

II/^IIhsII'?!!/,!! where || ■ denotes the norm.+ Therefore, there exists a /t > such 
that 

(F.24) 



(dcf/0*(gradAl[][f/,,(C)]) 



< K 



for all C G K(j. Then, since 



C{t) 



Im ( Tr I [R^Ut, (C) - f/^ (tf ; C)R] (t; C)fiU{t; C)]) (F.25) 



4«nAS3[^,,(C)] V 

at a local extremum of J'}^\ all local extrema of J'}^^ must lie within the hypersphere 
< K,/a. From this perspective, the role of the shape function can be interpreted 

\{dcu,,nsu)\\^^ 



Here, (dct/i 



cut! 



denotes the operator norm [fi4]: || [dcUti] 



sup 

5C/eTU(Wc) 



\SU\\m 
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as moving undesirably-shaped controls outside of this hypersphere, to the region where 
they can no longer act as local extrema of J'}^^. Note that both the space K^^ and 
the control objective J'}^^ depend on the shape function ri{t). When a shape function 
is used in an iterative optimization routine via (F.23), this results in the minimization 
of a different (though conceptually related) objective on a different space than when 

v{t) = 1. 
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